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Abstract: The quantum marginal problem asks what local spectra are consistent with 
a given spectrum of a joint state of a composite quantum system. This setting, also re- 
ferred to as the question of the compatibility of local spectra, has several applications 
in quantum information theory. Here, we introduce the analogue of this statement for 
Gaussian states for any number of modes, and solve it in generality, for pure and mixed 
states, both concerning necessary and sufficient conditions. Formally, our result can be 
viewed as an analogue of the Sing-Thompson Theorem (respectively Horn's Lemma), 
characterizing the relationship between main diagonal elements and singular values of 
a complex matrix: We find necessary and sufficient conditions for vectors (di, . . . , d n ) 
and (ci, . . . , c„) to be the symplectic eigenvalues and symplectic main diagonal el- 
ements of a strictly positive real matrix, respectively. More physically speaking, this 
result determines what local temperatures or entropies are consistent with a pure or 
mixed Gaussian state of several modes. We find that this result implies a solution to the 
problem of sharing of entanglement in pure Gaussian states and allows for estimating 
the global entropy of non-Gaussian states based on local measurements. Implications 
to the actual preparation of multi-mode continuous-variable entangled states are dis- 
cussed. We compare the findings with the marginal problem for qubits, the solution of 
which for pure states has a strikingly similar and in fact simple form. 

1. Introduction 

What reduced states are compatible with a quantum state of a composite system? The 
study of this question has in fact a long tradition - as the natural quantum analogue 
of the marginal problem in classical probability theory. Very recently, this problem, 
now coined the quantum marginal problem, has seen a revival of interest, motivated 
by applications in the context of quantum information theory jlHEElUMl. In fact, in 
the quantum information setting, notably in quantum channel capacity expressions, in 
assessments of quantum communication protocols, or in the separability problem, one 
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often encounters questions of compatibility of reductions with global quantum states 



Since it is only natural to look at the full orbit under local unitary operations, the 
quantum marginal problem immediately translates to a question of the compatibility of 
spectra of quantum states. The mixed quantum marginal problem then amounts to the 
following question: Is there a state p of a quantum system with n subsystems, each with 
a reduction pk, that is consistent with 



for k = l,...,n,r and denoting the respective vectors of spectra. In the pure 
marginal problem, one assumes p = \ip)(ip\ to be pure. In the condensed-matter con- 
text 11211131 . related questions are also of interest: For example, once one had classified 
all possible two-qubit reductions of translationally invariant quantum states, then one 
would be able to obtain the ground state energy of any nearest-neighbor Hamiltonian 
of a spin chain. The quantum marginal problem was solved in several steps: Higuchi et 
al. U solved the pure quantum marginal problem for qubits. Subsequently, Bravyi was 
able to solve the mixed state case for two qubits, followed by Franz J6) and Higuchi 
[[2] for a three qutrit system. The general solution of the quantum marginal problem 
for finite-dimensional systems was found in the celebrated work of Klyachko (0, see 
also Refs. 01411151 . This is indeed a closed-form solution. Yet the number of constraints 
grows extremely rapidly with the system size, rendering the explicit check whether the 
conditions are satisfied unfeasible even for relatively small systems. 

In this work, we introduce the Gaussian version of the quantum marginal problem. 
Gaussian states play a key role in a number of contexts, specifically whenever bosonic 
modes and quadratic Hamiltonians become relevant, which are ubiquitous in quantum 
optical systems, free fields, and condensed matter lattice systems. For general infinite 
dimensional systems the marginals problem may well be intractable. However, given 
that in turn these Gaussian states can be described by merely their first and second 
moments H16II17I . one could reasonably hope that it could be possible to give a full 
account of the Gaussian quantum marginal problem. This gives rise, naturally, not to a 
condition to spectra of quantum states, but to symplectic spectra, as explained below. 
For the specific case of three modes, the result is known ||T8l . see also Ref. lH9l . In this 
work we will show that this program of characterizing the reductions of Gaussian states 
can be achieved in generality, even concerning both necessary and sufficient conditions. 
This means that one can give a complete answer to what reductions entangled Gaussian 
states can possible haveQ 

Equivalently, we can describe this Gaussian marginal problem as a problem of com- 
patibility of temperatures of standard harmonic systems: Given a state p, what local 
temperatures - or equivalently for single modes, what local entropies - are compatible 
with this joint state? Of course, one can always take the temperatures to be equal. But 
if they are different, they constrain each other in a fairly subtle way, as we will see. In 
a sense, the result gives rise to the interesting situation that by looking at local tem- 
peratures, one can assess whether these reductions may possibly originate from a joint 



We refer here to the marginal problem for Gaussian states, which are quantum states fully defined by their 
first and second moments of canonical coordinates. However, clearly, our result equally applies to general and 
hence non-Gaussian states, in that it fully answers the question what local second moments are consistent with 
global second moments of quantum states of several modes. 



f7ll8ll9lfl0l[rai. 



spec( ( o) = r, 
spec(/9 fe ) = r fe 



(1) 

(2) 
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Fig. 1. Solution of the Gaussian marginal problem. The set of possible reductions with symplectic spectra 
(ci, . . . , c n ) of a correlated or entangled Gaussian state p with symplectic spectrum (di, . . . , d n ) is char- 
acterized by the given remarkably simple necessary and sufficient set of n + 1 inequalities. For example, 
from local measurements, one can hence infer about the consistency with the purity of the joint state. It also 
governs the sharing of correlations in Gaussian states. 



system in a pure state. Finally, it is important to note, since sufficiency of the condi- 
tions is always proven by an explicit construction, the result also implies a recipe for 
preparing multi-mode continuous-variable entangled states. 



2. Main result 

We consider states on n modes, and consider reductions to single modes. Gaussian 
states are represented by the matrix of second moments, the 2n x 2n covariance matrix 
7 of the system, together with the vector /i of first moments. For a definition and a 
survey of properties, see Refs. II16II17I . In this language, the vacuum state of a standard 
oscillator becomes 7 = I2, as the 2x2 identity matrix. The canonical commutation 
relations are embodied in the symplectic matrix 



fc=i 



1 

-1 



(3) 



for n modes. The covariance matrices of n modes are exactly those real matrices satis- 
fying 

7 + ia > 0, (4) 

which is simply a statement of the Heisenberg uncertainty principle. The first moments 
can always be made zero locally, and are hence not interesting for our purposes here. 
Note also that the set of Gaussian states is closed under reductions, so reduced states of 
Gaussian states are always Gaussian as well. 

Real matrices that leave the symplectic form invariant, SaS T = a, form the real 
symplectic group Sp(2n, R). In the same way as symmetric matrices M can be diago- 
nalized with orthogonal matrices to a diagonal matrix OMO T = D, one can diagonal- 
ize strictly positive matrices using such S £ Sp(2n, R), according to 

SMS T = D. (5) 

The simply counted main diagonal elements of D form then the symplectic spectrum 
of M, and the collection of symplectic eigenvalues can be abbreviated as sspec(TU) = 
(dx,.. .,dn), 

D = diag(di,di, . . . ,d n ,d n ). (6) 
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This procedure is nothing but the familiar normal mode decomposition. In turn, by 
definition, the symplectic eigenvalues are given by the square roots of the eigenvalues 
of the matrix —Ma Ma. Again, for the vacuum, the symplectic eigenvalues are all 
given by unity. In a mild abuse of notation, we will refer to the symplectic spectrum of 
a Gaussian state as the symplectic spectrum of the respective covariance matrix. 

Finally, for a given covariance matrix 7, and in fact any strictly positive real matrix, 
we refer to the symplectic main diagonal elements (ci, . . . , c„) as the symplectic eigen- 
values of the 2x2 main diagonal blocks. This is the natural analogue of main diagonal 
elements. Equivalently, the symplectic main diagonal elements are the main diagonal 
elements after the main diagonal 2x2 blocks have been brought into the form 



7k 



c k 
c k 



(7) 



We are now in the position to state our main result, see Fig. 1. It relates the symplectic 
spectrum of composite systems to the ones of the reductions. We will first state it as a 
mere matrix constraint, then as the actual Gaussian marginal problem, and finally for 
the important special case of having a pure joint state. 

Theorem 1 (Necessary and sufficient conditions). Let (di,...,d n ) and (ci, . . . , c„) 
be two vectors of positive numbers in non-decreasing order. Then there exists a strictly 
positive real 2n x 2n-matrix 7 such that (di , . . . ,d n ) are its symplectic eigenvalues and 
(ci, . . . , c n ) the symplectic main diagonal elements if and only if the n + 1 conditions 

k k 

>^2dj, k = l,...,n (8) 

j=i 3=1 

n— 1 n—1 

On - ^ Cj < d n - ^2 dj (9) 

are satisfied. 

This set of inequalities may be conceived as a general analogue of the Sing-Thompson 
theorem [20,21 ,22), see below. More physically speaking, this means the following: 

Corollary 1 (Gaussian marginal problem). Assume that p is a Gaussian state of n 
modes satisfying sspec(p) — (di , . . . , d n ). Then the possible reduced states p k to each 
of the individual modes k = 1, . . . , n are exactly those Gaussian states with 



sspec(p k ) = c k (10) 



satisfying Eq. (O and (O. 



These conditions hence fully characterize the possible reduced marginal states. For 
two modes, n = 2, for example, the given conditions read 

c 1 + c 2 >d 1 +d 2 , (11) 
C2 — ci < tfe — di, (12) 

for C2 > C\ and d 2 > d\. The constraint ci > d\ is then automatically satisfied. 
For pure Gaussian states, the above conditions take a specifically simple form. Quite 
strikingly, we will see that the resulting conditions very much resemble the situation of 
the marginal problem for qubits. 
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Corollary 2 (Pure Gaussian marginal problem). Let p = \ip) (ip\ be a pure Gaussian 
state ofn modes. Then the set (b\ + 1, ...,&„ + 1) of symplectic eigenvalues 

sspec(pk) = b k + 1 (13) 

k = 1, . . . , n, of the reduced states p k of each of the n modes is given by the set defined 
by 

bj<J2 b k (14) 

for all j, for bj > 0. 

To reiterate, these conditions are necessary and sufficient for the local symplectic spec- 
tra being consistent with the global state being a pure Gaussian state. 

Equivalently, this can be put as follows: If 7 is the covariance matrix of a pure 
Gaussian state with reductions 



Ik 



b k + l 
b k + l 



(15) 



k = 1, . . . ,n. Then, Eq. (fT~4-b defines the local temperatures T k per mode consistent 
with the whole system being in a pure Gaussian state, according to 

b k = 2(exp(l/T fe ) - (16) 

for the standard harmonic oscillator (an oscillator with unit mass and frequency). The 
above condition hence determines the temperatures that modes can have, given that 
a composite system is in a pure Gaussian state. The form of Eq. (fl3T l can always be 
achieved by means of local rotations and squeezings in phase space. One can hence 
equally think in terms of local symplectic spectra or local temperatures. 

It is instructive to compare the results for the pure Gaussian marginal problem with 
the one for qubits as solved in Ref. HI. There, it has been found that for a system 
consisting of n qubits, one has 

A,<^A fc (17) 

for the spectral values r k = (X kl 1 — A^.), A^ G [0, 1]. Moreover, these conditions are 
both necessary and sufficient. It is remarkable that this form is identical with the result 
for n single modes 

h<J2h, (18) 

b k > 0, as necessary and sufficient conditions. Again, the admissible symplectic eigen- 
values are defined by a cone the base of which is formed by a simplex. Note that the 
methods used in Ref. JT] to arrive at the above result are entirely different. Once again, 
a striking formal similarity between the case of qubit systems and Gaussian states is 
encountered. 

Finally, from the perspective of matrix analysis, the above result can be seen as a 
general analogue of the Sing-Thompson Theorem 120112 1 Il22l (or Horn's Lemma 11231 in 
case of Hermitian matrices), first posed in Ref. 1241 . where the role of singular values 
is taken by the symplectic eigenvalues. 
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Sing-Thompson Theorem (120,21,221) Let (x\, . . . ,x n ) be complex numbers such 
that \xk \ are non-increasingly ordered and let (2/1, . . . , y n ) be non-increasingly ordered 
positive numbers. Then annx n matrix exists with X\ , . . . , x n as its main diagonal and 
2/1, . . . ,y n as its singular values if and only if 

k h 

k = l,...,n, (19) 

3=1 J=l 
n— 1 n— 1 

\ X 3 \~\Xn\ <J2V] - yn- (20) 

J=l i=i 

It is interesting to see that - although the symplectic group Sp(2n, R) is not a compact 
group - there is so much formal similarity concerning the implications on main diagonal 
elements of matrices. Note, however, that the ordering of singular values and symplectic 
eigenvalues, respectively, is different in Theorem 1 and in the Sing-Thompson theorem. 



3. Proof 

As a preparation of the proof, we will identify a simple set of necessary conditions that 
constrains the possible reductions that are consistent with the assumption that the state 
is pure and Gaussian. These simple conditions derive from a connection between the 
symplectic trace and the trace of the covariance matrix. Quite surprisingly, we will see 
that they already define the full set of possible marginals consistent with a Gaussian 
state of n modes. We shall start by stating the condition to the reductions. 

Lemma 1 (Symplectic trace). Let 7 be a strictly positive real 2n x In-matrix such that 
its main diagonal 2x2 blocks are given by Eq. HT^for c^ G [1, 00). Then the symplectic 
eigenvalues (di, . . . ,d n ) of the matrix 7 satisfy 

n n 

J2 d k<J2 Ck - (21) 

fc=i fe=i 

Proof: Note that the right hand side of Eq. (I2TI 1 is nothing but half the trace of the 
covariance matrix 7, whereas the left hand side is the symplectic trace str(7) of 7, so 

n 

str( 7 ) = 5^ (22) 
i=i 

if sspec(7) = (di, . . . , d n ), see, e.g., Ref. Il25l . We arrive at this relationship by making 
use of a property of the trace-norm. The symplectic eigenvalues d%, . . . , d n of 7 are 
given by the square roots of the simply counted eigenvalues of the matrix (za)7(«cr)7 
II16II17I . Hence, the symplectic spectrum is just given by the spectrum of the matrix 

M=| 7 1 / 2 ( 1( t)7 1 /2| j (23) 
where | • | denotes the matrix absolute value ll26l . So we have that 



n 

2^d fc =tr(A/) = || 7 1/2 ^)7 1/2 ||i, 
fc=i 



(24) 
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where || • ||i is the trace norm. The property we wish to prove then immediately follows 
from the fact that the trace-norm is a unitarily invariant norm: this implies that 

n 

2]T4 = |l7 1/ V)7 1/2 ||i< IIMtIIi, (25) 
fc=i 

as || AS || i < ||£> j4|| i for any matrices A, B for which AB is Hermitian. This inequality 
holds for any unitarily invariant norm whenever AB is a normal operator ll27l . Now, 
since any covariance matrix is positive, 7 > 0, and the largest singular value of ia is 
clearly given by unity, we can finally conclude that 

n n 

2^<4<|| 7 || 1 =tr(7) = 2£> fe) (26) 

fc=l k=l 

which is the statement that we intended to show. ■ 



This observation implies as an immediate consequence a necessary condition for 
the possible reductions, given a Gaussian state of an n-mode system: Let 7 be the 
covariance matrix of a Gaussian pure state of n modes, with reductions as above. We 
can think of the state as a bi-partite state between a distinguished mode labeled k, 
without loss of generality being the last mode k = n, and the rest of the system. We 
can in fact Schmidt decompose this pure state with respect to this split using Gaussian 
unitary operations [30 29 35 1. This means that we can find symplectic transformations 
S A e Sp(2(n - 1), R) and S B G Sp(2, R) such that 



(S A ®S B )7(S A ®S B ) q 



A C 
C T B 



(27) 



A = diag(l, . . . , l,a n ,a n ), 
B = diag (a n ,a„), 



where 

(28) 
(29) 

with some 2(n — 1) x 2-matrix C. The symplectic eigenvalues of modes 1, . . . , n — 1 
are hence given by 1, . . . , 1, a n . The above statement therefore implies the inequality 

n - 2 + a n < ai + ■ ■ ■ + a„_i, (30) 

or, by substituting bk = a,k — 1 for all k = 1, 



. n, 

b„ < h + b 2 + • • • + b n - 



(31) 



This must obviously hold for all distinguished modes and not only the last one, and 
hence, we arrive at the following simple necessary conditions: 

Corollary 3 (Necessary conditions for pure states). Let 7 be the covariance matrix of 
a pure Gaussian state with thermal reductions 



Ik 



, n. Then, for all j, 



bk + 1 
b k + l 



(32) 



(33) 
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That is, the largest value of bj cannot exceed the sum of all the other ones. 

So far, we have assumed the global state p to be a pure state. In the full problem, 
however, we may of course allow p to be any Gaussian state, and hence a mixed one, 
with symplectic spectrum 

sspec(p) = {d u ...,d n ) > (34) 

instead of being (1, . . . , 1). This is the Gaussian analogue of the mixed marginal prob- 
lem. For this mixed state case, we provide necessary conditions for the main reductions, 
in form of n inequalities on partial sums, and one where the largest symplectic eigen- 
value of a reduction plays an important role. The first set of n conditions is up to the 
different ordering a weak majorization relation for symplectic eigenvalues, which is in 
fact essentially a corollary of a result from Ref. Il28l due to Hiroshima. The second state- 
ment, the n + 1-th condition, as well as showing sufficiency of the general conditions, 
will turn out to be significantly more involved. 

Lemma 2 (Necessity of the first n conditions). Let (di, . . . , d n ), and (ci, . . . , c n ) be 
defined as in Theorem[l] For any given (di, . . . , d n ), the admissible (ci, . . . , c„) satisfy 

k k 

for all k — 1, . . . , ri. 

Proof: Let S G Sp(2n, R) be the matrix from the symplectic group that brings 7 
into diagonal form, so 

SjS T = diag(di,di, . . . ,d n ,d n ). (36) 

The main diagonal elements of 7, in turn, again without loss of generality in non- 
decreasing order, are given by (cx, . . . , c n ). Now according to Ref. l28l , we have that 

k 

min tr(T7T T ) = 2^rfj (37) 
3=1 

for k = 1, . . . , n, where the minimum is taken over all real 2k x 2n-matrices T for 
which 

Ta n T T = a k . (38) 

Here, ak denotes the symplectic matrix on k modes as defined in Eq. (|3). Now we can 
actually take S G Sp(2n, R) according to S = 1, we see that T, consisting of the first 
2k rows of S, satisfies Eq. d38l . Since this submatrix does not necessarily correspond 
to a minimum in Eq. d37b . we find 

k k 

2^c J =tr(T 7 T T )>2^, (39) 

3=1 3=1 

for any k = 1, . . . , n. ■ 



We will now prove the necessity of the n + 1-th inequality constraint in TheoremQ] 
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Lemma 3 (Necessity of the last condition). Let [d\,...,d n ), and (ci , . . . , c„) be de- 
fined as in Theorem^ For any given vector of symplectic eigenvalues [d\, . . . , d n ), the 
admissible (ci, . . . , c n ) satisfy 

n—l n— 1 

Cn-^cj <d n ~J2 d r ( 4 °) 

i=i j=i 

Proof: We will define the function / : <S„ — > R, where 5„ is the set of strictly 
positive real 2n x 2n-matrices, as follows: We define the vector c = (a, . . . , c„) as 



j 



(722-1,23-1723,23 -72j-l,2j) 1/2 ( 41 ) 



j = 1 , . . . , n, as the usual vector of symplectic spectra of each of the n modes, and then 
set 

n 

/( 7 ):=2max(c)-^c j . (42) 

j'=i 

For a diagonal matrix D = diag(d l7 d\, . . . , d n , d n ) with entries in non-decreasing 
order, we have 

n-l 

f(D) = d„-J2dj. (43) 

i=i 

We will now investigate the orbit of this function / under the symplectic group, 

/ = sup {x G R : x = f(SDS T ), S E Sp(2n, R)} , (44) 

and will see that the supremum is actually attained as a maximum for 5=1. Each of 
the matrices 7 = SDS T have by construction the same symplectic spectrum as D. This 
is a variation over 2n 2 + n real parameters, as any 5 G Sp(2rt, R) can be decomposed 
according to the Eider decomposition as 

5 = OQV, (45) 

where 0,V E K(n) := Sp(2n, R) n 0(2n) and 

Q G {(«!, 1/zi, ...,Z„, 1/Zn) : z k G R\{0}} . (46) 

That is, O, V reflect passive operations, whereas Q stands for a squeezing operation. 

We will now see that the maximum of this function / - which exists, albeit the 
group being non-compact - is actually attained when the matrix is already diagonal. 
This means that in general, we have that 

/ = 2 maxsspec(7) — str(7). (47) 

For any global maximum, any local variation will not increase this function further. 
Let us start from some 7 = SDS T . For any such co variance matrix 7 we can find a 
T G Sp(2(n - 1),R) such that 



(T © l 2 h(T © I2? 



E F 
F T G 



--■ Y, (48) 
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where 



E = diag(ci, ci, . . . , d n _ x , c' n _ x ) 



(49) 



is a (2n — 2) x (2n — 2) matrix and G is a 2 x 2 matrix. Using Lemma[T| again, we find 
that 

71—1 n—1 



so 



/(V) > /(7)- 



(51) 



In other words, it does not restrict generality to assume the final covariance matrix to 
be of the form as in the right hand side of Eq. ((48), and we will use the notation 



7 = SDS 1 = 



E F 
F T G 



(52) 



withS = diag(ci,ci, . . . ,c„_i,c„_i) and G = diag(c„,c„). 

We can now investigate submatrices of 7 associated with modes m and n, 1 < m < 

n, 

Cm I2 Cn,r 



M„ 



(53) 



This we can always bring to a diagonal form, using symplectic diagonalization, only 
affecting the main diagonal elements of modes n and m, and leaving the other main 
diagonal elements invariant. This brings this submatrix into the form 



Mi 



c' m l 2 
c' n l 2 



with c' n > c' m . From Lemma [5] we know that 

d — c' < r 



(54) 



(55) 



so we have increased the function /, whenever C n , m 7^ 0. Hence, for global and hence 
local optimality, we have to have C„. m = 0. However, each of the matrices C„, m = 
for all m = 1, . . . , n — 1 exactly if the matrix 7 is already diagonal. 

What remains to be shown is that the function / is bounded from above, to exclude 
the case that the maximum does not even exist. One way to show this is to make use 
of the upper bound in Lemma|4]to have for every covariance matrix 7 with symplectic 
spectrum (di, . . . , d n ) in non-decreasing order 



/(7) <E^' + ( 3 - 2n ) d i 



(56) 



j=2 



which shows that / is always bounded from above. If 7 is merely a strictly positive 
real matrix, but no covariance matrix, an upper bound follows from a rescaling with a 
positive number. ■ 



We now prove the upper bound required for the proof of Lemma[3] 
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Lemma 4 (Upper bound). Let (di, . . . , d n ), and (ci, . . . , c n ) be defined as in Theorem 
\l\ and 7 be additionally a 2n x 2n covariance matrix. For any given (d\ , . . . , d n ), the 
admissible (ci, . . . , c„) satisfy 



n — 1 n 

^c, < Vrfj ■ + (3 - 2n)d 1 , 



/ , ^3 
3 = 1 



/ , 3 
3=2 



Proof: We start from a An x 4rt-covariance matrix 



7 



A C 
C T A 



corresponding to a pure Gaussian state, where 

n 



dfe 
d k 



fc=i 



o 



-(4 - 1) 1/2 



(57) 



(58) 



(59) 



(60) 



are real 2n x 2n-matrices. Physically, this means that we start from a collection of n 
two mode squeezed states, with the property that the reduction to the first n modes is 
just a diagonal covariance matrix with symplectic eigenvalues {d-y, . . . , d n ), again in 
non-decreasing order. Let us first assume that di = 1, this assumption will be relaxed 
later. Let us now consider 



Si 
1 



1 



ST 
1 



SiASf SiC 
C T S? A 



(61) 



for Si G Sp(2n, M). Obviously, the set we seek to characterize is the set B of main 
diagonals of the upper left block 

U = SiASf (62) 

of this matrix. We can always start from a diagonal matrix having the symplectic eigen- 
values on the main diagonal, and consider the orbit under all symplectic transformations 
S G Sp(4n,M). 

We will now relax the problem by allowing all S G Sp(4n, R) instead of symplectic 
transformations of the form S = Si © 1, Si G Sp(2n, R). We hence consider the 
full orbit under all symplectic transformations. This set C D B is characterized by the 
reductions to single modes of 



i = S-/S 7 



A' C 
C' T A 

-- diag(di,di, . 



(63) 

,dn.,d n ) This includes 



for some S G Sp(4n, R), such that again A 
the case (|6TT >. 

We are now in the position to make use of the statement that we have established 
before: From exploiting the Schmidt decomposition on the level of second moments, 
and using Lemma Q] relating the trace to the symplectic trace, we find 



n — 1 n 

°3 ^ J2 d J 



3 - 2n, 



(64) 
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as d\ — 1 was assumed. 

Let us now consider the case of di > 1. We will apply the previous result, after 
appropriately rescaling the covariance matrix. Indeed, we can construct a covariance 
matrix 7 as in Eq. ( 1581 for 

(di, . . . ,d n ) = (l,d 2 /di, ■ ■ • ,d„/di). (65) 

We then investigate the orbit of 7 under the symplectic group, and look at the main 
diagonal elements of SjS 7 '. By construction, we have that 7 + ia > 0. We can hence 
apply Eq. (l64t to this case. Multiplying both sides of Eq. (l63l l by d\ gives rise to the 
condition in Eq. (l57l i. ■ 



Lemma 5 (Solution to two-mode problem). There exists a strictly positive real 4x4- 
matrix 7 with main diagonal blocks diag{c\ , C\ ) , diag(c 2 , c 2 ) and symplectic eigenval- 
ues (di, c? 2 ) if and only if 



c-i + c 2 > di + d 2 , 
C2 — ci < d 2 — di, 



(66) 
(67) 



assuming c 2 > Ci a«d d 2 > c?i- Moreover, c\ — c 2 = c?i — c? 2 if and only the 2x2 o/f 
diagonal block ofj vanishes. 

Proof: The necessary conditions that |ci — c 2 | < |c?i — c£ 2 are a consequence of 
Lemma H The necessary conditions c\ + c 2 > c?i + <i 2 and ci > c?i have been previ- 
ously shown in Lemma|2] Hence, we have to show that these conditions can in fact be 
achieved. This can be done by considering a 



ci e 

Oci / 

e c 2 

/ c 2 



Sdmg(d 1 ,d 1 ,d 2 ,d 2 )S 1 



(68) 



The relationship between ci, c 2 , e, / and di, d 2 is given by 

d? /2 = (c? + C 2 + 2e/ 

± (cf + 4 + Aefcj - 2c\{c\ - 2ef) + 4 Cl c 2 (e 2 + / 2 )) 1/2 )/2, (69) 



as di , d 2 are the square roots of the eigenvalues of — (770-7 fTBll . compare also Ref. 
An elementary analysis shows that the above inequalities can always be achieved. Also, 
the extremal values are achieved if and only if e = / = 0. ■ 

What we finally need to show is that the conditions that we have derived are in fact 
sufficient. This will be the most involved statement. 



Lemma 6 (Sufficiency of the conditions). For any vectors (ci, . . . , c„) and (di , . . . , d n ) 

satisfying the conditions (0 and (0 there exists a 2n x 2n strictly positive real matrix 
with diag(c\, . . . , c n ) as its symplectic main diagonal elements and (di, . 
symplectic eigenvalues. 



,d n ) as its 
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Proof: The argument will essentially be an argument by induction, in several ways 
resembling the argument put forth in Refs. [20 .21 22]. The underlying idea of the proof 
is essentially as follows: On using the given constraints, one constructs an appropriate 
2(n — 1) x 2(n — 1) -matrix, in a way that it can be combined to the desired 2n x 2n- 
matrix by means of an appropriate S £ Sp(4, R) acting on a 4 x 4 submatrix only. 
Note, however, that compared to Ref. 1221 . we look at variations over the non-compact 
symplectic group Sp(2n, R), and not the compact U(2n). 

For a single mode, n = 1, there is nothing to be shown. For two modes, Lemma 
[5] provides the sufficiency of the conditions. Let us hence assume that we are given 
vectors (ci, . . . , c„) and (d%, . . . , d n ) as above, and that we have already shown that 
for 2(n — 1) x 2(n — l)-matrices, the conditions (|8]i and © are indeed sufficient. We 
complete the proof by explicitly constructing an 2n x 2n-matrix with the stated property. 

We have that c\ > d\ by assumption. We could also have c\ > dj for some 2 < j < 
n, so let k £ {1, . . . , n} be the largest index such that 

ci > d k . (70) 

Let us first consider the case that k < n — 2, and we will consider the cases k = n — 1 
and k = n later. Then we can set x := dk + d k+ i — cy, which means that x > 0, and 
that all conditions 

ci +x > dk + d k+ i, (71) 
ci - x > d k - dk+i, (72) 
-ci + x > d k - d k +i (73) 

are satisfied: (T7TT > by definition, ( 1721 because c\ > d k and ( 1731 as d k +i > C\. This 
means that we can find a matrix of the form 



7 



cil 2 C 
C T xl 2 



(74) 



for some 2 x 2-matrix C, with symplectic eigenvalues (dk,dk+i), using Lemma [5] 
Therefore, the matrix 

7" = 7' ediag(di,di,d 2 ,rf2, ■ ■ ■ , <4-i, d k+2 , d k +2, ■ ■ - ,d n ,d n ) (75) 

has the symplectic spectrum (d±, . . . , d n )- 

We will now show that we can construct a 2(n — 1) X 2(n — 1) matrix 7"' with 
symplectic eigenvalues (o?i, . . . , d k _i,x, d k + 2l ■ ■ ■ , d n ) and main diagonal elements 
(C2, C2, . . . , c„, c„), by invoking the induction assumption. This matrix 7"' we can in- 
deed construct, as we have 

C2 + --- + C; > rfi H hrfj-i, / = 2,...,fc, (76) 

C2 H hc fc+ i>(iiH 1-4-1+2:, (77) 

c 2 + h c s > di -\ h d k -i + x + d k+2 H \- d s ,s = k + 2, . . . , nl{78) 

as one can show using d k < c\ < d k+ i and x = d k + d k+ i — c\. Also, we have 

c n - c 2 c„_i <d n -di d k —i x (4+2 d n , (79) 

fulfilling all of the condition that we need invoking the induction assumption to con- 
struct 7"'. This matrix has the same symplectic eigenvalues as the right lower 2(n — 
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1) x 2(n — 1) submatrix 7"" of 7". Therefore, there exists an S G Sp(2(n — 1), R) 
such that 

7"" = Si" S T . (80) 

So the matrix 

7:= (1 2 ©S) 7 "(1 2 ®S) T (81) 

has the symplectic spectrum (di, . . . , d n ) and symplectic main diagonal elements (ci, . . . 
Hence, by invoking the induction assumption, we have been able to construct the de- 
sired matrix with the appropriate symplectic spectrum and main diagonal elements. 
Note that only two-mode operations have been needed in order to achieve this goal. 

We now turn to the two remaining cases, k = n and k = n — 1. In both cases 
this means that we have c\ > d n -i, as d n > d n -\, and both cases can be treated in 
actually exactly the same manner. Obviously, this implies that also c n > c\ > d n -i- 
We can now define again an x, by means of a set of inequalities. This construction is 
very similar to the one in Ref. l22l . We can require on the one hand 

x > max{d„_i, rf„_i + d n - c n ,d n -i - d n + c„, 

di H h d n - 2 + Cn^! - ci c„_ 2 }- (82) 

On the other hand, we can require 

x < min{d n - d„_i + c„, c\ H V c„_i - d± d«-2, 0}. (83) 

Both these conditions can be simultaneously satisfied, making use of c„ > c„_i and 
c n > dn-i. This in turn means that we have 

c n + x > d n -i + d n , (84) 
c n -x> d„_i - d n , (85) 

X - C„, > d n -l - d n , (86) 

where the latter two inequalities mean that \x — c n \ < |c? n -i — d n \. Moreover, we satisfy 
all the inequalities 

cH +Q>cH \-di,l = l,...,n-2, (87) 

ci + h c n _i > di H h d„-2 + x, (88) 

and 

C n -1 — C\ C„_ 2 < X - d\ d n -2- (89) 

Again, we can hence invoke the induction assumption, and construct in the same way as 
before the desired covariance matrix with symplectic spectrum (di, . . . ,d n ) and sym- 
plectic main diagonal elements (c\, . . . , c n ). This ends the proof of sufficiency of the 
given conditions. ■ 
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4. Physical implications of the result and outlook 

The results found in this work can also be read as a full specification of what multi- 
partite Gaussian states may be prepared: Since the argument is constructive it readily 
provides a recipe of how to construct multi-mode Gaussian entangled states with all 
possible local entropies: For pure states, starting from squeezed modes, all is needed 
is a network of passive operations. Applied to optical systems of several modes, no- 
tably, this gives rise to a protocol to prepare multi-mode pure-state entangled light of 
all possible entanglement structures from squeezed light, using passive linear optical 
networks, via 

7 = OPO T , (90) 

with P = (zi, 1/zi, . . . , z n , l/z n ), Zk G R\{0}, and O G K(n). P is the covariance 
matrix of squeezed single modes, whereas O represents the passive optical network. 
The latter can readily be broken down to a network of beam splitters and phase shifters, 
according to Ref. Il34l . Hence, our result also generalizes the preparation of Ref. IT8l 
from the case of three modes to any number of modes. Similarly, for mixed states, the 
given result readily defines a preparation procedure, but now using also squeezers in 
general. 

The above statement also settles the question of the sharing of entanglement of sin- 
gle modes versus the rest of the system in a multi-mode system: For a pure Gaussian 
state with d\ = ■ ■ ■ = d„ = 1, the entanglement entropy £yi{x n}\{j} °f a m °de 
labeled j with respect to the rest of the system is given by 

where s : [1, oo) — > [0, oo) is a monotone increasing, concave function. 

Corollary 4 (Entanglement sharing in pure Gaussian states). For pure Gaussian 
states, the set of all possible entanglement values of a single mode with respect to the 
system is given by 

{ E l\{2,. ..,«}, ■ • • , #n|{l,...,n-l}) G < { s i c l)i ■■, s ( c n)) ■ Cj - 1 < ^(Cfc - 1), Cj > 1 

(92) 

This result is an immediate consequence of the above pure marginal problem, Corol- 
lary In fact, this is for pure Gaussian states more than a monogamy inequality: it 
constitutes a full characterization of the complete set of consistent degrees of entangle- 
ment. 

A further practically useful application of our result is the following: It tells us how 
pure a state must have been, based on the information available from measuring local 
properties like local photon numbers. This is expected to be a very desirable tool in an 
experimental context: In optical systems, such measurements are readily available with 
homodyne or photon counting measurements. 

Corollary 5 (Locally measuring global purity in non-Gaussian states). Let us as- 
sume that one has o.cquired knowledge about the local symplectic eigenvalues Ci , . . . , c n 
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of a global state p. Then one can infer that the global von-Neumann entropy S(p) of p 
satisfies 

S{p)<s{y j cX (93) 

This estimate is true regardless whether the state p is a Gaussian state or not. 

Proof: Let us denote with to the Gaussian state with the same covariance matrix 
7 > as the (unknown) state p. The vectors (di, . . . ,d n ) and (ci, . . . , c n ) are the 
symplectic eigenvalues and symplectic main diagonal elements of u>, respectively. From 
the fact that diag(di , d\ , . . . , d n , d n ) reflects a tensor product of Gaussian states, we can 
conclude that 

n 

5M=^ S (d,). (94) 

i=i 

In turn, from LemmaQ]we find that 

n n 

5>>X>. (95) 

i=i i=i 



By means of an extremality property of the von-Neumann entropy (see, e.g., Ref. 
36 1) that a Gaussian state has the largest von-Neumann entropy for fixed second mo- 
ments, we find that S(p) < S(uj). Since the function s : [1, oo) — > [0, oo) defined in 
Eq. (l9Tb is concave and monotone increasing, we have that 

n 

S(p) < S K") ^ S ( d l + • • • + d n) < S(C1 + • • • + Cn). (96) 
3=1 

This is the statement that we intended to prove. Clearly, this bound is tight, as is 
obvious when applying the inequality to the Gaussian state with covariance matrix 
diag(di ,d 1 ,...,d n ,d n ) itself. ■ 

For example, if obtains c\ — 3/2 = C2 = 3/2 and C3 = 2 in local measurements 
on the local photon number, then one finds that the global state necessarily satisfies 
S(p) < s(5). This is a powerful tool when local measurements in optical systems are 
more accessible than global ones, for example, when no phase reference is available, or 
bringing modes together is a difficult task. 

To finally turn to the role of Gaussian operations in this work: Our result highlights 
an observation that has been encountered already a number of times in the literature: 
That global Gaussian operations applied to many modes at once are often hardly more 
powerful than when applied to pairs of modes. This resembles to some extent the situ- 
ation in the distillation of entangled Gaussian states by means of Gaussian operations 

HEEBIES]. 

In this work, we have given a complete characterization of reductions of pure or 
mixed Gaussian states. In this way, we have also given a general picture of the possi- 
bility of sharing quantum correlations in a continuous-variable setting. Since our proof 
is constructive, it also gives rise to a general recipe to generate multi-mode entangled 
states with all possible reductions. Formally, we established a connection to a com- 
patibility argument of symplectic spectra, by means of new matrix inequalities fully 
characterizing the set in question. These matrix inequalities formally resembles the 
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well-known Sing-Thompson Theorem relating singular values to main diagonal ele- 
ments. It is the hope that this work can provide a significant insight into the achievable 
correlations in composite quantum systems of many modes. 
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